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Abstract 

The Lorentzian Hamiltonian constraint is solved for isotropic loop 
quantum cosmology coupled to a massless scalar field. As in the Eu- 
clidean case, the discreteness of quantum geometry removes the clas- 
sical singularity from the quantum Friedmann models. In spite of the 
absence of the classical singularity, a modified DeWitt initial condi- 
tion is incompatible with a late-time smooth behavior. Further, the 
smooth behavior is recovered only for positive or negatives times but 
not both. An important feature, which is shared with the Euclidean 
case, is a minimal initial energy of the order of the Planck energy re- 
quired for the system to evolve dynamically. By forming wave packets 
of the matter field an explicit evolution in terms of an internal time 
is obtained. 

1 Introduction 

The recent paper reports on a quantization of a Euclidean cosmology 
in which a massless scalar field was coupled to a spatially flat Friedmann- 
Robertson- Walker model. Set in the framework of isotropic loop quantum 
gravity (ILQC) j2] , this model describes cosmological evolution in terms of a 
discrete time. This discreteness is a direct consequence of the kinematics of 
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loop quantum gravity which predicts that geometric quantities such as area 
lUEllH], volume lU El El El CD] , and angle [TT] have a discrete spectrum. In 
[2j the Hamiltonian constraint of quantum gravity, expressed in a discrete 
basis of volume eigenfunctions, and the matter Hamiltonian act on the state 
function of the coupled system. This model avoids the classical singularity. 
Due to the quantization of the inverse scale factor of the Friedmann model 
[TB] the total Hamiltonian constraint vanishes at the classical singularity. 
At small values of the volume the wave function displays a distinct discrete 
behavior. Nevertheless as the volume of the model grows, the wave func- 
tion approaches a continuous function which is a solution of the asymptotic 
Wheeler - DeWitt differential equation. The wave function thereby meets 
the semi-classicality requirements of quantum cosmology [12]. 

In the present paper we investigate the Lorentzian Hamiltonian constraint 
for the same massless scalar field model which is constructed along the pre- 
scriptions of the general theory [21], as it was introduced in 0. In the 
spatially fiat Friedmann model, however, the classical full Lorentzian con- 
straint is proportional to the Euclidean one, so there is a quantum ambiguity 
and one can, as in recent work [2^ , alternatively consider the Euclidean con- 
straint operator proportional to the full one. In the semiclassical regime the 
results of both versions converge. The choice in this paper is closer to the 
full theory of quantum geometry |21j . 

We concentrate on three aspects of the ILQC framework. First, the 
Hamiltonian constraint at early times leads to a consistency relation for ini- 
tial data - the "dynamical initial conditions" of ^j. Although the Lorentzian 
constraint is higher order than the Euclidean constraint, we find the same re- 
lation in the model. Second, to select an essentially unique solution one may 
require late time solutions to be smooth. In such a "pre-classical" [20] state 
the wave function at late times does not vary strongly on short intervals. We 
find that this criteria again selects an essentially unique solution. However, 
in contrast with the Euclidean model evolution backwards to negative time 
destroys the pre-classical condition. Third, an unexpected feature and, at 
the same time the main result of is the occurrence of a threshold for dy- 
namical evolution of the model. For the wave function to have a dynamical 
interpretation with respect to an internal time there is a minimal energy on 
the order of the Planck energy, concentrated initially in a volume of the size 
of a Planck volume. 

Scalar field quantum cosmology was first considered by Blyth and Isham 
[2]. Canonically quantizing the reduced model, they explored the dynamics 



2 



for several choices of time. The current work is in sharp contrast to this 
older work in that the discreteness of ILQG manifests itself in early times 
and completely changes the status of initial conditions of the cosmological 
model. 

As we are dealing with a generalization of the framework of PP in section 
12] we present the prerequisites very briefly . In section |21 the discrete Hamil- 
tonian constraint is solved numerically and the solutions are discussed. In 
section HI the asymptotic Wheeler - DeWitt equation and its solutions are 
constructed and compared with the Euclidean case. In section El an essen- 
tially unique wave function with a sufficiently smooth late time behavior to 
represent a classical universe, is singled out from the general solution of the 
Hamiltonian constraint. Section |B1 contains a study of possible wave packets 
in the scalar field and their classical interpretation. 



2 Prelimaries: Loop Quantum Cosmology 

We start with a quantization of the spatially fiat Friedmann-Robertson- 
Walker model with a scalar field as matter source. The metric 

represents a homogenous, isotropic model with one dynamical degree of free- 
dom, represented, for example, by the scale factor a{t) and its canonically 
conjugate momentum. In the framework of loop quantum gravity, applied to 
isotropic models ^2], the metric variable a{t) is replaced by a triad variable 
p. This triad variable can assume both signs according to the two possible 
orientations of a triad. By defining 

p := sgn(a) (1) 

the domain of a now extends to negative values. The conjugate momentum 
is 

c = ^d. (2) 

Isotropic loop quantum cosmology yields a discrete basis \n) of volume 
eigenstates, 

3 

V\n) = (^jll) ' ^(\n\-l)\n\{\n\ + l) \n) =: Vi(|„|_,)|n), (3) 
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where the integers n label discrete values of p (or the scale factor) 

6 

lip 



n = -^p. (4) 



In this basis the inverse scale factor operator, which is constructed inde- 
pendently from a, is also diagonal |12j, 



H^|n) = 16(7/^)-2 (yvT^- ^Vi|,|_i) \n). (5) 

This, in contrast to the inverse of the scale factor operator, is a densely 
defined operator on the Hilbert space spanned by the basis {|^)}. 

For the purposes of evaluating the Hamiltonian constraint the most im- 
portant feature of this construction is that the eigenvalue zero of the volume 
is threefold degenerate, as can be seen from Q, it vanishes in |0) and in 

I ± 1), whereas the eigenvalue of is zero only on |0). The latter state 
assumes the role analogous to the classical singularity a = 0. The vanishing 
of the inverse scale factor is a pure quantum feature, in sharp contrast to the 
divergence of the classical a"^ for a = 0. 

An arbitrary state \s) of the cosmological model can be expressed as 

oo 

\s) = ^Sn\n). (6) 

— oo 

In the ILQC framework the full Lorentzian Hamiltonian constraint H\s) = 
assumes the form of a difference equation for the coefficients s„ j2] 



''n+8'^n+4'^n+8 



3(7«^p)"^ i(l + 7"^) sgn(n + 8) (Vi|„+8| - '^i|n+8|-i) K+A. 

- sgn(n + 4) (Vi |„+4| - Vi |„+4|_i) s„+4 

-2sgn(n) {Vi\^\ - Vi|„|_i) (|(1 + 7~2)(A;-A;++4 + k+k~_^) - 

- sgn(n - 4) (Vi |„_4| - |„_4|_i) s„_4 (7) 

+ 7"^) sgn(n - 8) (Vi|„_8| - Vi|„_8|_i) k~_8k-_^Sn^8 = 0. 

Thanks to the sign factors the coefficient sq of the singular state drops out 
and can be set equal to zero. In this sense the model based on loop quantum 
gravity is singularity- free. 
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y\{\n+l\-l) - Vi(|„_3|_i) j ( Vi|„_3| - Vi|„_3|_i + Vi|„+i| - 1^i|„,+i| 



The fc's, coming from the extrinsic curvature contribution to the "kinetic" 
term of the Hamiltonian, have the following expressions in terms of volume 
eigenvalues 

K = mir' 

and 

K = 3(7^ 



-1 



Vi 



i(|n+5|-l) - l^i(|n+l|-l) 



^i|n+l| - '^i|n+l|- 



1 + V, 



Vi 



2(|n+3|-l) ^i(ln-lhl) 



Vi|„__i|_i + Vi|„_^3| - \/i|„+3|_i 



i|„+5| - Vi\n+5\- 
Vi 



and are subject to the identities 



t+ — h~ 



and 



t+ — Ic~ 



(10) 
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These identities make only one kind of /c's necessary; we choose k~^. 

The model under consideration includes a massless scalar field 0. Its 
classically Hamiltonian in the Friedmann- Robert son- Walker metric is Hj^ = 



^p'ia ^. When quantized it becomes P 
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(12) 



which is the energy operator per unit coordinate volume. Here the quantiza- 
tion of |a|~^ (0) is crucial for it renders the classical Hamiltonian finite when 
a = 0. In the total Hamiltonian couples with a sign factor sgn(n) to the 
gravitational part. 

As in |1J we assume the states of the coupled system are of the form 



(9) 



\s) = ^Sn{(j)) |n), 



(13) 



where the dependence on (f) is contained in the coefficients s„ of the quan- 
tum geometry basis vectors \n). In our calculations we assume that the 
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0-dependence of the state vector is given by an eigenfunction x of the matter 
Hamiltonian, characterized by cj, 

Sn(0)=:s„x.(0):=s„e*^<^. (14) 

As we see in the next section, this ansatz resuhs in a finite difference equation 
for the model. 

3 The Hamiltonian constraint equation 

The difference equation resulting from the Hamiltonian constraint ((7j) is of 
order 16 so that the naive expectation is that the solution should contain 
16 free parameters. The solution decomposes into four independent series of 
coefficients, namely 54^+1, i = 0, ... ,3, m G Z, with four free parameters 
each. In the series with i = Sq drops out, because its coefficient in the 
gravitational hamiltonian contains sgn(O) and H(f,{0) is the zero operator by 

virtue of its construction in terms of \a\~^ So, for example, for n = 8, 

the Hamiltonian constraint ((Tj) with the matter Hamiltonian, does not relate 
So to S4, . . . , S16, but instead gives a consistency condition for S4,...,Si6, 
reducing the number of free parameters from four to three. This i = series 
is considered as fundamental. By applying the "pre-classicality condition" 
at late times we are able to pick out essentially unique solutions to the other 
series. The idea is that, with the fundamental series, the other three series are 
selected by smooth interpolation. Thus, the 16 free parameters are reduced 
to 3. 

To facilitate the handling of the difference equation we simplify it by 
introducing rescaled coefficients 

tn := (Fih - 5„ (15) 

and use the abbreviation 

/3:=i(l + 7-')- (16) 

To avoid confusion, we remind the reader that Sn denote the full coefficients, 
Sn the 0-independent parts and t„ the rescaled Sn used for solving the Hamil- 
tonian constraint 

{H + sgn(n)i^^) ^ sM\s) = 0. (17) 
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The sign factor in front of the matter Hamiltonian as well as the signs in (|7j) 
give rise to a relative sign between the two (classically disjoint) sectors with 
p > and p < (see p]). They render the above constraint equation time 
symmetric, provided Htf){n) = H^{—n). For the moment, we consider only 
positive values of n. 

By inserting n = 8 into the total Hamiltonian constraint equation [H + 
H^) s„((/))|s) = we obtain tiQ in terms of t^, ts and ti2, 



tl6 



12 



2048cj2 (v/T^ 



2-f3[{k 



127 



(18) 



ts + ^4 



(This is the above-mentioned consistency condition.) The general form of 
the difference equation for n > 12 is 



'-n+4 



n+4 



2 - /3 m^,? + iK?] + 



(19) 



tn + tn_4 l3k^_^k^ ■ tfi-g 



Equation (fT^ is complicated enough (the worst complications are hidden 
in the fc's) so that numerical methods are useful. In figure Q "we plot the 
numerical solution for the initial conditions = 1, tg = tu = and a 
parameter u such that /t-y^/iu = 10^- ^^is numerical solution we use the 

value of 7 = ~ 0.13 for the Immirzi parameter.^ Although the solutions 
were obtained from the initial values ^4 = 1, tg = '^12 = the solution is 
generic in that the solution is qualitatively the same when the initial values 
are varied by as much as 10^. 

In contrast to the Euclidean case, persistent short-wavelength oscilla- 
tions continue to late times - even deep into what we would expect to be 



"'^This value merges from isolated horizons work of 18 and in a general model of grav- 
itational statistical mechanics |19j . 
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the semiclassical regime. These short-wavelength oscillations appear to be 
superimposed on a smooth function of the same shape as the Euclidean wave 
function. 




Figure 1: A solution Sm of (fnijl with initial values t4 = 1, = ti2 = 0, 
7 = 0.13, and m = 4n so that < n < 800. The points Sm of the wave 
function are connected by lines to clearly display the oscillatory character. 



4 The continuum limit 

The physical implications of the Lorentzian Hamiltonian constraint can be 
deduced without a detailed solution of (fTIJ|) . but from the asymptotic smooth 
mean-value functions, which are approximated for large values of n (or the 
volume). These functions are solutions of a Wheeler-DeWitt differential 
equation, which can be derived from as continuum limit. To distin- 
guish between discrete and continuous quantities we return the continuous 
variable p, whose relation to n is given by (^. 
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We also need some asymptotic expansions: 



8 96 
- ~ 24-i(7/|)tn^ = 17/1^, (21) 

^-^n^)'- (Ti8)M7/|)^--i = (h^|)%-i. (22) 

Assuming the 0-dependence to be given by the function Xui.'P) of equation 
(|T^. we rewrite the Hamiltonian constraint (fT^ in the form 

(^^n+8^n+4 '^n+S ^ {^t+i) +(^n) ''^n + ^"-8 j (23) 

6 



- (W4 - 2t„ + t„_4) - -^^^ Vi„ - Vi„_i 
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so that we have two kinds of differences on the left-hand side, differences of 
products of t's by /c's and differences of Vs alone. In the continuous limit for 
large n the later ones are second derivatives with respect to p. Thanks to 
the fact that the fc's are, to order equal to one, the former differences 

are also approximated by the second derivative in leading order. Hence the 
continuous limit of the left-hand side becomes 

in which we write the asymptotic form of t„ as t{p). (For the moment we 
pretend that n is continuous.) Together with the leading term of the right- 
hand side this gives a Cauchy-Euler equation of order 2 

t(p) =ple±^^i'^P, (25) 



with the solutions 
where VL is given by 



Vt = \J^Kun-p^ - I. (26) 
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From this solution we may construct the continuous hmit S;^{p,(f)) of 
from the relation s„(0) = tn/ (Ki/2 — Vn/2-1) Xui{(p)- Taking into account that 
Vi„ — V^n-i S^^s as in leading order we have 

s^(p,0) = e±*^'"^e±^t^. (27) 

The continuous limit of the wave function is thus the same as in the Euclidean 
case The same is true for the existence of the critical value ujcrit = 
The behavior of the wave function is crucially different according to the value 
of u. The wave function only displays the asymptotic oscillatory behavior 
when u > cjcrit- The value cUcrit corresponds to a threshold value of an initial 
energy of the order of the Planck energy, concentrated in a volume of the order 
Ip at the classical singularity (see P).^ Only if such an amount of energy 
or more is present, dynamical evolution of the model is possible. Figure |21 
shows the approximation of an asymptotic wave function of the type (j27|) by 
a solution of the Hamiltonian constraint (|T^ for = 0. 



5 Pre-classicality 

5.1 Finding the pre-classical solution 

In j2ni it is argued that in loop quantum cosmology there always exists a 
"pre-classical" solution which is smooth at late times. To find this solution 
among the three-parameter family of the general solution, it is convenient to 
look at a different asymptotic limit than achieved in equation (|2l|l . 

This limit may be found in [SHj. But to keep our discussion self-contained 
we outline it here. The approximation is based on replacing the fc's by their 
asymptotic value of 1 and by considering the matter contribution to the con- 
straint equation as small and constant. Indeed, {^jVni2 — ^Jyufi^xf^Ynj^ — 
14/2-1 )~^ goes as for n ^ 1, so this approximation is valid on finite 
ranges of ra, determined by the condition that does not change much on 
such a range. We define 

_ 2048^^ (vA^-v/n^)° 

■^The precise value may be modified by quantum ambiguities |16[I17 |. 
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Figure 2: The numerical result (points) and the asymptotic solution (solid 
line) for the 0-independent part s^rn of the wave function. The conditions 
are the same as in figure ^ but here 100 < m < 1000. 

On this scale the high frequency oscillations are visible in the "cloud" of 
points around the asymptotic solution. 



which is much smaller than 1 for large n and use the index m = n/4 in the 
following. 

Considering P as constant we obtain an asymptotic difference equation 
with constant coefficients instead of a differential equation, 

/?W2 - Wi + 2(1 - /? + P)tm - t^^i + (3t^-2 = 0, (29) 

conserving discreteness. An exponential ansatz tm = e^"^^ leads to a quadratic 
equation for cos-i? with the approximate solutions 

P 1 — 2(3 P 
cos^9o~l + Y3^ and cost^i ^ T^^' ^^^'> 

From this we obtain 4 solutions, e^^"^^° and e^*™''^^ , where 'do goes asymptoti- 
cally to zero and "i^i does not. In consequence, e^*™'^° approximate a constant 
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in a range of m which is much smaller than the wavelength of the main os- 
cillation. On an arbitrary long range these solutions approximate a smooth 
solution of (j23)- This solution is called pre-classical. This analysis shows that 
the asymptotic behavior of the three linearly independent functions obtained 
from the three initial parameters is of the form 

+ a» cos dim + h^^ sin dim, i = 1,2, 3, (31) 

where /wdw ^ Teal, p (respectively n)-dependent part of a solution of the 
asymptotic Wheeler-DeWitt equation. For any linearly independent triple of 
such functions there is obviously a unique, nontrivial linear combination in 
which the sine and cosine of dim cancel. 

In the numerical solution extreme fine tuning of the initial values is nec- 
essary to find smooth solutions.^ It is better to determine these values indi- 
rectly. With our set of 3 initial values we can produce a variety of solutions, 
whose average curves, as shown in figure 2, are in arbitrary phase relations 
with one another. These smooth functions form the 2-dimensional space of 
solutions of the asymptotic differential equation (j2H). 

To determine the pre-classical solution we evolve smooth initial data for 
sufficiently large values of n backwards to n = 0. For this purpose we may 
choose any set of 4 smooth data; every set gives rise to a solution of (j24p in 
the asymptotic region. Back evolution yields fairly smooth functions back 
to a early domain in n. In this domain and earlier the functions begin to 
oscillate rapidly and, in general, they do not fulfill the consistency condition 
(fTSj) for ^4 to ti6 and thus for sq to Siq. 

Constructing two linearly independent asymptotic solutions, we obtain 
a basis of solutions of the asymptotic Wheeler-DeWitt differential equation, 
from which it is possible to construct that linear combination, which fulfills 
the consistency condition. This combination yields the initial data §4, ss, and 
S12, whose forward evolution, shown in figure 3, is the pre-classical solution 
for the given values of 7 and the field energy. 

As a practical matter the essential uniqueness of the pre-classical solution 
has to be understood within a certain latitude. Late time solutions are very 
sensitive to small changes in the initial data for n = 4, 8, 12. For instance, 
the solution depicted in figure 3 is not the optimal one, although it appears 
sufficiently smooth in the figure. The solution arises as a linear combination 

■^Changing the initial values by one part in 10^^ will re-introduce the oscillatory "quan- 
tum foam" at later times. 
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Figure 3: The pre-classical solution for the same values of 7 and the field 
energy as in figure 1 for n from 40 to 2000. For smaller values of n the 
amplitude increases, the largest coefficient S4 is equal to 150222.0. 

of two pre-classical functions, evolved back from the late time initial values 
(si988,Si992,Si996,S2ooo) equal to (1,1,1,1) and (0.7,0.8,0.9,1), respectively. 
In this domain of n these linearly combined initial data are a sufficiently 
good approximation to pre-classical wave functions so that the resulting 
graphic looks perfectly smooth. By a slight modification, they could be 
adjusted to a smooth asymptotic solution of the form in order to improve 
the pre-classicality of the resulting function. 

One might hope that a symmetry principle might select the pre-classical 
wave function. For instance we might consider data which is symmetric or 
antisymmetric about the time of the classical singularity. As we see in the 
next section though, neither of these cases select a pre-classical solution. 

5.2 Evolution through the classical singularity 

The absence of singularities in loop quantum cosmology enables the wave 
function to evolve through the time of the classical singularity into a domain 
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with negative values of n. Evolving the pre-classical wave function backwards 
reveals a further essential difference between the Euclidean and Lorentzian 
cases: When n approaches zero from the positive side, the amplitude of the 
wave function begins to increase. Beyond zero it continues growing and after 
large and rapid oscillations the wave function settles down to the superpo- 
sition of the two oscillations, already familiar from figures ^ and |21 Note, 
however, that for negative n's the wavefunction has a very large amplitude 
compared to positive n's. For our values of 7 and the field energy, the ratio 
between the amplitudes for negative and positive values of n is of the order 
of 10^^. Naively, the wave function for negative and positive n suggests that 
there is vanishing probability of the pre-classical cosmology at positive n. 



-1000 



-800 



-600 



^00 



-200 



Figure 4: Continuation of the above pre-classical wave function to —1000 < 
m < —1. 

Near the classical singularity, Lorentzian wave functions display large os- 
cillations. The larger the energy of the scalar field, the longer is the "quan- 
tum regime" of these oscillations at positive n's and the later the pre-classical 
behavior sets in. 
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6 Wave packets and dynamics 



So far, the functional dependence of the wave function on the field was as- 
sumed to be of the form exp(±2^0). The field being massless and spatially 
constant, (j) is formally equivalent to a configuration variable of a particle 
with gravity acting as a one-dimensional potential. In this analogy the criti- 
cal energy separates free states from bound ones (although we did not define 
a measure on n to give the wave functions a proper meaning of a probability 
amplitude at a certain value of 0). 

To investigate wave packets in 0, made from superpositions of x^'s for 
different cj's, we use 



Sn{(t>) = sn- / dcue-^('^-'^°)'e-*^<^ (32) 



with ojo ^ LUcrit and A ^ (u;o — ^^crit)^^- The two conditions assure that the 
contribution from uj < ujcrit is negligible and no under critical wave functions 
are included. 

The Wheeler - DeWitt equation being linear, we obtain in the continuous 
limit a superposition of functions (j77j) 



s(p,0) = (27rAn2)-3 rdcue-^(--'^«)%-*^<^e^v(^)'-'"^^ 



oo 



With normalization in (f) with respect to the natural inner product 

/■oo 

{s\t)= / d0s*(p,0)t(p,0) 



this leads to the following modulated GauB function in 

1 1/1 11 p </'\2' /Irip <p\ 

s{p,(j)) ^ {27rXh^)—ie~^^^^t'^> e"^"^^^"*^ (33) 

with a maximum at = AlEE ^nd a width 2\/X. 

Considering now p as time variable, we can calculate p-dependent expec- 
tation values of (p and the (kinetic) field energy. For (p 

(0(p)) = 7r^lnp, (34) 



15 



showing a growth during the expansion of the universe. This is an indication 
that in a more sophisticated model with an inherent notion of particles, like a 
massive or a spatially non-constant field, one could expect particle creation. 
The energy expectation value 

E^p) := {H^ip)) = (-|p-t ^) = I {ul + 1/4A) (35) 

shows a decreasing field energy per unit coordinate volume. Remembering 
the relation between the discrete counterpart n of p and the scale factor 
a we find that the energy scales as V~^. The energy density per physical 
volume goes therefore as V~'^. 

In the above interpretation, when the radius of the universe is considered 
as internal time, there is no problem with superpositions of expanding and 
contraction universes or universes going backwards in time, as it would be, 
when acted as time. So the choice of n or p seems to be much more 
natural and less problematic - the wave functions and expectation values of 
the matter field evolve with the scale factor. 

7 Conclusions 

Comparing the present Lorentzian with the Euclidean Hamiltonian con- 
straint, one notices three similarities: (1) The asymptotic continuous wave 
function is of the same form. (2) There is the same minimal initial energy 
(up to quantum ambiguities in both cases ^Hl) for dynamical evolution of 
the wave function. (3) The extension of the quantum regime between the 
classical singularity and the beginning of a smooth, pre-classical evolution 
increases with growing matter energy. Therefore the semiclassical limits of 
the wave functions coincide. This fact confirms that the difference between 
the two versions of the constraint is a quantum ambiguity of the spatially 
flat Friedmann model. 

In the quantum regime, the models differ essentially from each other. The 
Euclidean wave function is small near the singularity, thus satisfying what 
one could call a "modified DeWitt initial condition". The original DeWitt 
condition ensures that the wave function in standard quantum cosmology 
stays away from the singularity. The Lorentzian wave function, on the other 
hand, although being equal to zero at the singular state, is large in the 
immediate vicinity of the classical singularity. 
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In both cases it is possible to extend the principle series beyond the sin- 
gularity to negative values of the discrete internal time n. In the Euclidean 
case there is only one solution (up to a normalizing factor), which is smooth 
and symmetric with respect to n. In the Lorentzian case the unique^ smooth 
solution is embedded in a three-parameter family of Planck-scale varying so- 
lutions. Furthermore, pre-classicality is only one-sided, i.e. the wave function 
is pre-classical, at best, for n > or n < 0. 

Thus, if we assume that evolution from arbitrary negative to arbitrary 
positive values of the internal time is possible, we come to the following sce- 
nario. Whenever a pre-classical universe is to emerge from the considered 
model, it is preceded by a non-pre-classical "fuzzy" universe. The latter one 
contracts to a state of zero volume and bounces off the singularity. For very 
special initial conditions the wave function is reflected almost completely 
and a very small fraction, which will behave pre-classically, enters the do- 
main of positive time. In other words, provided appropriately adjusted initial 
conditions, the singularity acts as a "filter" , which keeps back all "fuzzy" im- 
purities and transmits only the purely pre-classical wave function. From the 
large ratio of the amplitudes of the fuzzy and the smooth wave functions and 
from the fact that there is no pre-classical wave function at all unless the 
parameters are very fine tuned, it appears that a (pre-) classical universe is 
very unlikely to emerge from our model. However, our intuition on probabil- 
ities, transition and reflection does not have a reliable framework. To obtain 
a sound notion of probability and, perhaps, of unitary evolution, one needs 
a suitable inner product on the space of wave functions. 

The absence of singularities and the possibility to evolve wave functions 
through a state with zero volume formally solves a problem of classical cos- 
mology and "standard quantum cosmology", where the scale factor of the 
universe is considered as continuous. Nevertheless our model does not have 
a clear interpretation. In fact it raises much the same kind of questions as 
standard cosmology. These questions, however, are rephrased in the context 
of ILQC in a way which might prove productive. 

Above all, there is the question of the beginning of the universe. Is a (pre- 
) classical, expanding universe preceded by a time-reversed, shrinking one? 
To discuss this question, one has to think about whether in an emphatically 
non-pre-classical evolution the interpretation of the parameter n as internal 

*Up to a small uncertainty, see the remark at the end of section 5.1. This is in addition 
to the normalizing factor. 
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time makes sense. Does the internal time pass from large negative values of n 
towards zero, or should we "sec" the beginning of the universe at the classical 
singularity? Do the positive and negative branches of the wave function 
describe two possible evolutions of the universe away from the initial zero- 
volume state? Does the huge ratio between the amplitudes of the non-pre- 
classical and the pre-classical wave function indicate only a tiny probability 
for a pre-classical universe to be created? 

The present model poses some interesting questions for ILQC. One short- 
coming (shared with all homogeneous models) is that the massless, spatially 
constant scalar field cannot imply the notion of particles, so that it is insuffi- 
cient for the description of cosmological particle production. Another is the 
issue of stability. It would be interesting to investigate whether these models 
would be stable under the inhomogeneous perturbations. 
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